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A  procedure  is  presented  for  calculating  the  electromagnetic  field  in 
the  vicinity  of  an  infinitesimally  thin  conducting  circular  cylindrical  shell 
with  an  infinitely  long  slot  illuminated  by  a  TE  plane  wave.  This  field  is  the 
short-circuit  field  plus  the  field  due  to  the  slot.  The  short-circuit 
field  is  the  field  that  would  exist  if  there  were  no  slot.  The  field  due 
to  the  slot  stems  from  the  <j>  component  of  the  electric  field  in  the  slot. 
Approximated  by  a  linear  combination  of  expansion  functions  that  satisfy  the 
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JO.  ABSTRACT  (continued)  / 

^  edge  conditions,  this,  field  component  is  determined  by  testing  the  equation 
of  continuity  of  the/axial  magnetic  field  in  the  slot  with  non-negative 
powers  of ^  where  |^|  is  the  angular  distance  from  the  center  of  the  slot. 

Alternative  expansion  functions  are  introduced  to  preserve  the 
accuracy  of  the  solution  in  the  vicinity  of  an  internal  resonance. 
Containing  the  resonant  electric  field,  the  first  alternative  expan¬ 
sion  function  is  normalized  so  that  its  axial  magnetic  field  remains 
finite  as  the  resonant  frequency  is  approached.  The  resonant  electric 
field  is  that  interior  electric  field  whose/fx component  approaches 
zero  on  the  closed  surface  consisting  of  the  cylindrical  shell  and  the 
slot  as  the  frequency  approaches  the  resonant /frequency .  The  rest  of 
the  alternative  expansion  functions  are  orthogonal  to  the  resonant 
electric  field.  / 
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I.  INTRODUCTION 


Scattering  of  a  TE  electromagnetic  wave  by  a  conducting  circular 
cylinder  with  an  infinitely  long  axial  slot  was  computed  in  [1]  and  [2] 
by  three  different  methods:  numerical  solution  of  the  aperture  field 
Integral  equation,  numerical  solution  of  the  electric  field  integral  equa¬ 
tion,  and  numerical  solution  of  the  magnetic  field  integral  equation.  The 
method  of  numerical  solution  of  the  electric  field  integral  equation  was  also 
used  in  [3].  When  the  slot  is  so  narrow  that  the  field  inside  the  cylinder  is 
a  small  fraction  of  the  incident  field,  the  field  inside  the  cylinder  can  not 
be  accurately  determined  by  the  method  of  numerical  solution  of  the  elec¬ 
tric  field  integral  equation  because,  in  this  method,  the  field  inside  the 
cylinder  is  the  sum  of  the  incident  field  and  the  field  due  to  the  numeri¬ 
cally  calculated  electric  current  on  the  surface  of  the  cylinder.  The 
incident  field  is  the  field  that  would  exist  if  the  cylinder  were  absent. 
Similarly,  numerical  solution  of  the  magnetic  field  integral  equation  will 
not  accurately  determine  the  field  inside  the  cylinder  when  this  field  is 
a  small  fraction  of  the  incident  field.  Furthermore,  the  magnetic  field 
integral  equation  applies  only  to  a  conducting  cylindrical  shell  of  finite 
thickness.  The  scattering  problem  was  also  formulated  as  a  dual  series 
problem  and  solved  by  techniques  borrowed  from  the  Riemann-Hilbert  problem 
in  complex  variable  theory  £  4 ]  —  [8 ] . 

In  the  present  report,  we  solve  the  problem  of  scattering  by  an 
infinitesimally  thin  conducting  circular  cylindrical  shell  with  a  narrow 
but  infinitely  long  slot  by  means  of  the  aperture  field  integral  equation. 

As  in  [2],  we  use  the  method  of  moments  to  obtain  a  numerical  solution  of 
the  aperture  field  integral  equation.  In  contrast  to  the  pulse  expansion 
functions  and  point  testing  functions  used  in  [1],  we  use  expansion 


functions  that  satisfy  the  edge  conditions  and  test  with  non-negative 
powers  of  $  where  j<}>|  is  the  angular  distance  from  the  center  of  the 
slot  aperture.  In  comparison  with  [9,  pp.  1392-1398]  where  only  one 
expansion  function  is  used,  the  aperture  field  integral  equation  is 
point  matched,  and  no  numerical  calculations  are  performed,  we  use  up 
to  4  expansion  and  testing  functions  and  calculate  fields  at  the  center 
of  the  cylinder  and  at  the  center  of  the  slot  aperture. 

In  our  formulation,  the  elements  of  the  moment  matrix  are  evalu¬ 
ated  by  replacing  each  expansion  function  by  a  very  large  but  finite 
number  of  terms  in  its  Fourier  series  with  respect  to  Bessel  func¬ 
tions  of  large  order  are  approximated  by  means  of  Debye's  asymptotic 
expansions  [10, p.  366].  To  obtain  a  moment  matrix  that  remains  well- 
conditioned  as  the  frequency  approaches  a  resonant  frequency  of  the 
associated  cavity,  we  choose  expansion  functions  so  that  the  magnetic 
field  of  only  one  of  them  contains  the  resonant  magnetic  field.  The 
associated  cavity  is  the  cavity  enclosed  by  the  complete  circular  cylin¬ 
drical  surface.  The  complete  circular  cylindrical  surface  is  the  con¬ 
ducting  shell  with  the  slot  covered  by  a  conducting  surface  whose  radius 
of  curvature  is  the  same  as  that  of  the  shell.  The  magnetic  field  of  an 
aperture  electric  field  expansion  function  is  the  magnetic  field  that 
exists  inside  the  cylinder  when  the  tangential  electric  field  vanishes 
on  the  conducting  circular  cylindrical  shell  and  is  equal  to  the  expan¬ 
sion  function  in  the  slot.  The  resonant  magnetic  field  is  the  magnetic 
field  associated  with  the  resonant  electric  field.  The  resonant  electric 
field  is  the  interior  electric  field  whose  $  component  approaches  zero  on 
the  complete  cylindrical  surface  as  the  resonant  frequency  is  approached. 


So  constructed,  our  moment  solution  remains  accurate  when  the  associated 


y 


cavity  is  resonant.  A  computer  program  was  written  to  implement  our  solu¬ 
tion.  This  computer  program  will  be  described  and  listed  in  a  forthcoming 
report . 


II.  MATHEMATICAL  FORMULATION 

The  formulation  for  the  TE  excited  cylindrical  shell  is  similar  to 
that  for  the  TM  excited  cylindrical  shell  in  [11].  Furthermore,  the  com¬ 
puter  program  that  was  written  to  implement  it  is  similar  to  the  one 
in  [12].  An  infinitesimally  thin,  perfectly  conducting  circular  cylin¬ 
drical  shell  of  radius  a  having  a  slit  aperture  of  half  angle  4>q  is 
illuminated  by  a  TE  polarized  plane  wave.  In  terms  of  the  polar  coordi¬ 
nates  (p,4>,z) ,  the  equation  of  the  shell  is  (p«*a,  4>o  <_  4>£  2tt— d)Q)  and  the 
incident  magnetic  field  H*  is  taken  to  be 


u  JkPcos(<Hx) 


— ze 


(1) 


where  a  time  factor  e',Wt  has  been  assumed  and  suppressed.  See  Fig.  1. 

The  incident  magnetic  field  is  the  magnetic  field  that  would  exist  if 

the  shell  were  absent.  In  (1),  u  is  the  unit  vector  in  the  z  direction, 

— z 

a  is  the  angle  of  the  direction  from  which  the  incident  wave  comes,  and 
k  ■  oo/ye  is  the  wave  number  where  y  and  e  are,  respectively,  the  per¬ 
meability  and  permittivity  of  the  medium  surrounding  the  shell.  Our  ob¬ 
jective  is  to  find  the  electric  and  magnetic  fields  everywhere  in  Fig.  1. 

We  close  the  aperture  with  an  infinitely  thin  curved  conducting 
strip  whose  equation  is  (p**a,  |4»|  £  $  )  and, as  shown  in  Fig.  2,  we  place 
the  magnetic  current  M  on  the  right-hand  side  of  the  closed  aperture  and 
-M  on  the  left-hand  side  of  the  closed  aperture.  These  magnetic  currents 
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L.  The  geometry. 
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2.  Equivalent  Situation. 


restore  the  tangential  electric  field  in  the  aperture.  The  situation  in 
Fig.  2  will  be  equivalent  to  that  in  Fig.  1  if  the  tangential  magnetic 
field  is  continuous  across  the  closed  aperture.  Thus,  since  the  magnetic 
field  of  a  TE  wave  is  in  the  z  direction,  we  require  that 


IT(-M)  -  Ht(M)  +  H®C  , 
z  —  z  —  z 


P  =  a 


1*1  1  *r 


(2) 


where  u  H  (-M)  is  the  magnetic  field  due  to  -M  placed  on  the  region  b 
z  z  — 

side  of  the  closed  aperture  in  Fig.  2,  u^  (M)  is  the  magnetic  field  due 
to  M  placed  on  the  region  a  side  of  the  closed  aperture  in  Fig.  2,  and 

SC 

u  H  is  the  magnetic  field  that  would  exist  if  the  incident  wave  impinged 
z  z 

on  the  conducting  shell  with  its  aperture  closed,  i.e.,  on  the  conducting 
circular  cylinder  of  radius  a.  Region  a  is  the  region  for  which  p  >_  a. 
Region  b  is  the  region  for  which  p  <_  a. 

Equation  (2)  is  recast  as 


Hb(-M)  -  Ha(M)  =  HSC  , 
z  —  z  —  z 


P  *  a 


1*1  i  *D 


(3) 


We  let 


M  =■  u  (M®  +  M°) 


(4) 


0  o 

where  M  is  even  in  4>  and  M  is  odd  in  <f>.  Because  the  magnetic  field 

0  Q 

due  to  u  M  is  even  in  $  and  that  due  to  u  M  is  odd  in  d>,  (3) 

— z  — z 

decomposes  into 


Hb(-u  M6)  -  Ha(u  Me)  =  H®Ce, 
z  — z  z  — z  z 


p  =  a 


*  I  1*0 


(5) 


TT^  f  wO  \  ,.3  /  »  .O  »  -SCO 

H  (-u  M)-H(uM)=H  \ 
z  — z  z  ■ — z  7  z 


1*0 


(6) 


where  H  is  the  even  part  of  H  and  H  is  the  odd  part  of  H 
z  z  z  z 


Seeking  an  approximate  numerical  solution  to  (5),  we  let 
4 


where 


M6  -  l  V®M( 
j-1  J  ' 


e  e 


(7) 


|<t>|  <  <P 


(8) 


otherwise 


The  expansion  functions  (8)  are  reasonable  because  the  equivalent  magnetic 

current  for  a  narrow  but  infinitely  long  slot  in  a  planar  conducting  screen 

/~2  2 

w  -  x 

[~ 2  o 

and  x/ J w  - 


x  terms  [13,  eq.  (44)]  where  w  is  one  half  of  the  slot  width 
and  |x|  is  the  distance  from  the  center  of  the  slot.  Substitution  of  (7) 
into  (5)  gives 


l  V®[Hb(-u  M*f)  -  Ha(u  Me)]  =  H1 
j.l  J  z  -z  j  z  -z  j 


see 

z 


p  =  a 


(9) 


|4>l  1  <t>0 


1)  6 

In  (9),  H  (— u  M  )  is  the  z  component  of  the  magnetic  field  in  region  a 

Z  Z  J 


d  0 

of  Fig.  3  and  H  (u  M. )  is  the  z  component  of  the  magnetic  field  in  region  b 
z  — z  j 


of  Fig.  3.  The  magnetic  fields  in  Fig.  3  are  obtained  by  solving  the  boundary 


value  problem  in  which  E^,  the  <}>  component  of  electric  field,  is  continuous 


at  p  =  a  where  it  is  given  by 


4>  j 


(10) 


In  terms  of  the  elementary  wave  functions  described  in  [14  ,  Sec.  5-1], 


Kn-mr, 


A-«h') 

Z  — Z  j 


»  S®  H^2^  (kp)  cos  (n<t>) 
1  v  Jn  n  _ 

n  n-0  H(2) ' (ka) 

n 


00  S®  J  (kp)  cos  (n<j>) 

I  V  J.n.  .  J1 _ _ _ _ _ 

n  n=0  J'(ka) 


where  n  =  /y/e  and  S®  is  an  unknown  coefficient. 

jn 

CL4  ,  eq.  (5-19)]  that 


It  is  evident  from 


E4>  “  jn  9  (kp) 


where  E.  is  the  d)  component  of  the  electric  field  in  Fig.  3  and  H  is  the 
q>  z 

z  component  of  the  magnetic  field  in  Fig.  3.  Substituting  (11)  into  (12) 
and  using  (10) ,  we  obtain 


M®  -  l  S®n  cos  (n0) 


From  (13),  S®  is  the  Fourier  coefficient  of  M®  given  by 
jn  3 


E  f 

Sjn  "  2^  j  cos  d<*> 


Here,  e  is  Neumann's  number 
n  r 


1  ,  n  «  0 


2  ,  n  >  1 


Substitution  of  (8)  into  (14)  gives 


S®  -  — 
jn  2tt 


(<t>/4>  )2^  2cos  (n$)  dcf> 


-4>0  yj  1  -  (4>/4>0>‘ 


The  coefficient  S.  is  evaluated  in  Appendix  A. 


W'i 

W 


rk 
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Using  (11)  and  the  Wronskian  [14,  eqs.  (D-12)  and  (D— 1 7 ) ] 


J  (ka)  H(2)'(ka)  -  J’(ka)  H(2)(ka)  *  - 
n  n  '  7  n  n  '  uka 


we  obtain 


v  „  00  S.  cos  (n<j>) 

[Hb(-u  M*)  -  Ha(u  M*)J  -  l  - IS - 7^y\l - 

2  "*  J  2  "2  J  0-a  -mka  niQ  Jl(ka)  hC2) 

n  n 


which,  when  substituted  into  (9),  gives 


_  4  00  S.  cos(n4>) 

_2_  y  ve  (  y  _ ^ _ ) 

*nka  j-1  i  n-0  J'(ka)  H(2)'(ka) 

n  n 


rHSce]  ,  | (|i|  <  <j> 

[  z  J  p=a  1 Y 1  -  ro 


It  is  evident  from  [14,  eq.  (5-116)]  that 


sc  ..jl  I  jVrcti).. 

1  a  Jo-a  irka  m'..  . 


2  P-.  vk.  n.LH(2)'(ka) 

n 


The  even  part  of  (20)  is  given  by 


°°  e  1  cos  (na)  cos  (n<£) 

rHscei  -  -  A L  y  _ _ 

1  *  "ka  nt0  2H(2>’(ka) 

n 


Substitution  of  (21)  into  (19)  and  subsequent  multiplication  by  Trnka/4 


00  S  cos  (n<t>) 


1  “T  w  .  V.WO  Vllt  / 

\  l  V?  (  I  - & - (-2yr- '  )--3n  I 

j-1  J  n-0  J'  (ka)  (ka)  n-0 


00  e  j  cos  (na)  cos  (n<J>) 
r  n 


2H(2) '(ka) 
n 


,  I  |£4>0 


Seeking  an  approximate  numerical  solution  to  (6), 


we  let 


M°  -  l  V°  M° 
j-1  J  3 


where 


m 


«*'**/ 

'm 

M 

fr'M 


m 

W 

$; 

m 


m 

m 


p 

Ess 


vVS 

vvX 

VV* 


m 


fv-  .  Jl  z’» 
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,  |  <p  |  <  4»o 


"I 

i 

,  otherwise 

Proceeding  as  In  the  development  (9)  -  (22),  we  obtain 


(24) 


.n 


,  4  °°  S  sin  (n<J>)  c  ........... 

r  I  - 72V - )  -  -  in  I  j.  »  ;M  (25) 

n-  '  n  (ka>  nml  Hn  <k*)  _ 

n  n  n 


0  L  i  t 

j«l  J  n-1  Jj^(ka)  H 
where 


jn  it 


-A 


(<(>/<) )  )^  *sin  (ncf>)  dt{) 


A7 


(26) 


(<t>/<PQ)' 


The  coefficient  S^n  is  evaluated  in  Appendix  A. 

When  the  frequency  is  very  close  to  a  resonant  frequency  of  the 
associated  cavity,  there  is  an  integer  p  such  that  |j^(ka)|  is  extremely 
small.  In  this  case,  it  can  be  seen  by  the  following  reasoning  that  any 
numerical  method  of  solution  will  fail  to  accurately  determine  .  The 
four  cos  (p<(»)  terms  on  the  left-hand  side  of  (22)  become  so  large  that 
none  of  the  other  terms  can  be  discerned.  As  a  result,  the  only  informa- 
that  can  be  extracted  from  (22)  is  that 


I  vT  s' 

j-i 


,e  ge 

j  Jp 


(27) 


This  single  equation  is  not  sufficient  to  determine  vf,  V®,  V®,  and  vf. 

12  3  4 

There  should  be  no  difficulty  in  solving  (25)  numerically  when 
|j^(ka)|  is  very  small  because  J^(ka)  does  not  appear  in  (25).  However, 


•  1 


11 


if  |jp(ka)|  is  very  small  where  p  is  a  positive  integer,  then  numerical 
solution  of  (25)  will  be  just  as  difficult  as  that  of  (22). 


III.  ALTERNATIVE  EXPANSION  FUNCTIONS 


In  this  section,  the  expansion  functions  are  changed  in  order  to 

eliminate  the  indeterminacy  encountered  in  the  last  two  paragraphs  of 

Section  II.  If  |j'(ka)|  is  extremely  small,  (7)  is  replaced  by 
P 


where 


M6-  l  V®M® 
j-1  2  J 


^  -  j;(ka) 

*j  “  +  ^  ^  ’  j-2’3’4 


where 


_e  .  e 
C  -  -  S.  /S, 
i  jp  ip 

Comparing  (28)  with  (7)  and  using  (29),  we  obtain 


V*-J'(ka)V*+  l  C JVJ 
w  i-2 


Ve  -  V® 

j  j 


j-2,3,4 


In  view  of  (30),  substitution  of  (31)  into  (22)  gives 

0  6 
-  S,  cos(p<jn  00  S,  cos  (n<f>) 

\  K  (*l2V - +  JP(ka)  2  - 5 — (2 y - > 

H'  '  (ka)  p  n-0  J'(ka)  IT  '  (ka) 

P  n*.  " 


I  «?<  I 

j-2  3  n-0  J' 

nj*p  n 


00  (S®  +  C®  S®  )  cos  (n<f) ) 

r  jn  j  In 


'(ka)  H(2)?(ka) 
n  n 


P 

Kw 
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„  e  j  cos  (na)  cos  (n<f>) 
00  n 


-  jn  l  9H(2), 


n-0  2H--0ca) 


.  I  <J>  I  £  4>q 


Numerical  solution  of  (32)  should  not  be  difficult  when  J*  (ka) I  is  very 

P 

small  because,  in  contrast  to  (22),  (32)  does  not  contain  any  huge  cos(p$) 
terms  which  obscure  all  the  other  terms. 

If  |j^(ka)|  is  very  small,  (25)  is  retained.  However,  if  |j^(ka)| 
is  very  small  where  p  is  a  positive  integer,  then  (23)  is  replaced  by 


where 


where 


M°=  I 

j=l  3  3 


%  - 


-  M°  +  C°  M°,  j-2,3,4 


C°  -  -  S°.  /S° 
j  JP  IP 


Comparing  (33)  with  (23)  and  using  (34),  we  obtain 


=  J'to,  V-  ♦  l  c°v‘ 


vj  -  •  J-2’3'4 


In  view  of  (35),  (36)  causes  (25)  to  be  replaced  by 


,  S?  sin(p4>)  00  S°  sin  (n$) 

1  r.  ('r  +  J'(ka)  I  -is - ^ - : 

2  1  H(2)  (ka)  P  n-1  J ' (ka)  H(2)  (ka) 

P  n*>  n 

i  4  00  (S°  +  C°S°  )  sin  (n4>) 

1  r  ?,o  ,  v  in  j  In'  x 


*j  l  ^  <  I  -J£ — 4? 

Z  j-2  3  n-1  J'  (ka) 

<_  n  n 

nj*p 


■-JH  i  J"31"  <"?>  3ln  (n*>  .  1*1  ;  * 

n-1  Hu;  (ka)  “  *■ 

n 


C, 'ij 'i 

IVl  iV 

Jvvi 

t«V 

WJi'i 

& 

?•:& 


| 

m 


■  > 

8 

ww 

IJ 

9 

« 

is 


i 

1 

Ij 

i 

M 


m 

to: 


Equation  (37)  is  valid  for  p  >  1.  If  p  ■  0,  we  revert  to  (25).  Numerical 


solution  of  (37)  should  not  be  difficult  when  |j^(ka)|  is  very  small  because, 
in  contrast  to  (25),  (37)  does  not  contain  any  huge  sin  (p<J>)  terms  which 
obscure  all  the  other  terms. 

Equations  (5)  and  (6)  can  be  viewed  as  operator  equations.  For 
example,  (5)  is 

L  M®  -  g  (38) 


where  L  is  an  operator  and  g  is  a  known  function.  Now,  (18)  is  LM^  where 
e 

Mj  is  given  by  (13).  Since  (8)  was  not  used  in  its  derivation,  (18)  can 
be  generalized  to  mean  that,  given  an  even  2tt  periodic  but  otherwise  arbi¬ 


trary  function  f(4>),  Lf  is  the  right-hand  side  of  (18)  where  {S®  }  are  the 


jnJ 


Fourier  coefficients  of  f.  Taking  f  -  cos  (p<j>) ,  we  obtain 


L(cos(p<f>) )  ■  - 2  c-°.3  .0*^  t -  (39) 

Ttnka  J'(ka)  IT  '  (ka) 

P  P 

When  ]jp(ka)|  is  very  small,  the  expansion  functions  {M®,  j *1 ,2,3,4}  are 
not  good  because  each  one  of  them  contains  cos  (p<p)  so  that  all  of  the 
functions  (L(M®),  j*l,2,3,4},  being  roughly  proportional  to  cos  (pd>) «  are 
nearly  indistinguishable  from  each  other.  The  expansion  functions 
(ft®,  j *1,2, 3, 4}  are  better  because  ft^,  which  contains  cos  (p4>) ,  is  normal¬ 
ized  so  that  Lft^  remains  finite,  and  none  of  the  functions  {ft®,  j*2,3,4} 

contains  cos  (p<j>) .  Therefore,  none  of  the  functions  {Lft®,  j*2,3,4}  con¬ 
tains  cos  (p<}>) .  Since  none  of  them  are  dominated  by  cos  (p<J>) ,  the 
functions  {Lft®,  j*2,3,4}  are  easily  distinguishable  from  each  other  and 
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f»  tt 


#.  ).•  a. 


IV.  TESTING  FUNCTIONS 


In  this  section,  4  linear  equations  are  extracted  from  the  functional 
equation  (32)  by  multiplying  (32)  by  each  of  4  even  functions  called  testing 
functions  and  integrating  from  to  with  respect  to  <f>.  Similarly,  four 
linear  equations  are  extracted  from  (37)  by  multiplying  (37)  by  each  of  four 
odd  testing  functions  and  integrating  over  the  slot. 

Called  (W®,  i*l,2,3,4},  the  even  testing  functions  are  defined  by 

W®  =  (4>/4>  )21"2  ,  i-1,2.3,4  (40) 

l  o 

0 

Multiplying  (32)  by  and  integrating  over  the  slot,  we  obtain  four  equa¬ 
tions,  one  for  each  value  of  i.  In  matrix  notation,  these  equations  are 

Ye$e  -  t*  (41) 

.g  0 

where  V  is  the  column  vector  jwhose  jth  element  is  and  Y  is  the 
4x4  matrix  whose  elements  are  given  by 


Fe  Se 
ip  ll 


e  e 

F7  si 
in  In 


Yn-  (2)>  p  +J^(ka>  l  - — 

(ka)  p  n-0  J'(ka)  IT  ;  (ka) 
P  n  n 


i=l ,2,3,4 


^  Fin(S1n  +  Ci  Sln> 


Ye  „  £  in  jn  j  in 

ij  n-0  J' (ka)  H(2)'(ka) 
,  n  n 

n?p 


i  -  1,2, 3, 4 


j  -  2,3,4 


where 


F®  -± 
in  2 


(<t>/4»  ) 2i_2cos (n$)d<f> 
o 


The  coefficient  F®n  is  evaluated  in  [11,  Appendix  B].  In  (41),  T®  is  the 


column  vector  whose  elements  are  given  by 


IJ--2JT,  l 


00  e  jnF®  cos  (na) 
v  nJ  in 
4  (j)  • 

n-0  2HV  }  (ka) 
n 


,  i-1,2.3,4 


sV  iVtv»v- 


I 

If 


i 


M 

p 

kt 

m 


I WiWS 


The  odd  testing  functions  are  defined  by 


W°  -  ($/<J>o)2i_1  ,  i  -  1,2, 3, 4 


(46) 


We  multiply  (37)  by  and  integrate  over  the  slot.  Doing  this  for  i-1,2,3, 
and  4,  we  obtain  the  matrix  equation 


o^o  -*o 
Y  V  -  1^ 


(47) 


*0  o 

where  V  is  the  column  vector  whose  jth  element  is  ,  and  Y  is  the  4x4 

matrix  whose  elements  are  given  by 


F°  S° 


Y,,  ■  +  l 


.O 

F  S 
in  In 


_H^a)  "  V~'  n-1  J’(ka)  H(2)’(ka) 
P  n 


,  i-1,2,3,4 


(48) 


•  F°  (S°  +  C°  S°  ) 

1°  m  ^  in  jn  j  In 

1J  n-1  J'(ka)  H(2),(ka) 


ni*p 


where 


in  2 


-<J> 


(♦/*  )21"1  sin  <„«  d* 


(49) 


(50) 


The  coefficient  F°n  is  evaluated  in  [11,  Appendix  B].  in  (47)  7°  is  the 


column  vector  whose  elements  are  given  by 


1°  “  "  2jn  l 


00  jn  F°  sin  (na) 
“  in 


n-1  H(2)'(ka) 

n 


i  -  1,2, 3, 4 


(51) 


Actually,  (47)— (49)  are  valid  only  for  p  >  0.  For  p  =  0,  we  multiply 
(25)  by  W°  and  integrate  over  the  slot.  Doing  this  for  i  =  1,2,3,  and  4, 
we  obtain  (47)  with  the  elements  of  ?°  given  by 


V°  -  V°  , 

J  J 


with  the  elements  of  Y  given  by 


(52) 


16 


Yij  -  l  — —  *  («> 

1J  n-l  J'(ka)  HU;  (ka) 
n  n 

and  with  the  elements  of  1°  given  by  (51). 

After  (41)  is  solved  for  $e,  Me  will  be  given  by  (7)  with 
(V®,  j*l,2,3,4}  given  by  (31).  After  (47)  is  solved  for  V°,  M°  will  be 
given  by  (23)  with  {V°,  J  “  1,2, 3, 4}  given  by  (36)  for  p  #  0  and, 
according  to  (52),  by 

V°  -  V°  ,  J  -  1,2, 3, 4  (54) 

0  O 

for  p  *  0.  Finally,  M  is  given  by  (4)  where  M  and  M  have  been  described 
in  the  previous  two  sentences: 


s  ’  *.  ^  (VH  +  W 


The  magnetic  field  in  Fig.  2  is  u  H  where 

— z  a 


Ha(M)  +  HSC  ,  p  >  a 
z  —  z 


Hbz(-M) 


,  p  <  a 


sc  a  b 

where  H  ,  H  (M) ,  and  H  (-M)  are  the  z  components  of  the  magnetic  fields 
z  z  z 

due  to  the  incident  wave,  M,  and  -M,  respectively.  From  (1)  and  [14, 
eqs.  (5-114)  and  (5-115)],  we  obtain 

oo  j » (ka) 

*  l  eJntJn(kc)  -  -TrT? -  H(2)  (kp)  ]  cos  (n(<D-00)  (f 

Z  n-0  n  n  H(2)  (ka)  n 

n 
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H  (M)  -  l  (V®  H  (u  MT)  +  V°H  (u  M?)) 
Z  -  j“i  j  z'-z  2  j  z  -z  2 


W.Hl.'t 


H°(-M)  -  l  (V®H*(-u  M*)  +  vV(-u  M°)) 


yz'^y  ]  v  vj 


where  H®(u  M®)  and  Hu(-u  M')  are  given  by  (11).  Similarly, 

Z  Z  J  Z  Z  J 


,b  /  »e 


00  s°  H(2)(kP)  sin  (nth) 


wos  j  v  jn  n  '  '  v 

VW  2,  ^ - Trp - 

2  ^  J  n  n=l  HU)  (ka) 

n 

_b,  o  i  c  S°„Jn(kp)  «ln  (b») 


so  that  (58)  expands  to 


a  1  *  „  ;  s!  «<2)<k0)  COS  (n*)  4  „  -  s”h'2) (kp>  sin  (n*) 

<(M)  -  i  (  I  v'(  I  ->-a  -T - )  +  I  v°<  I  Ja-a  - : 

n  j-1  J  n-0  H(2;  (ka)  j-1  J  n=l  H(2)  (ka) 


;*y, ('Jy 

M 


SM 


1 

i®: 


4  00  S  J  (kp)  cos  (n<{») 


H  (-M)  *  “T  [  I  V*(  l  -jn  n  j.fkaN - )  +  B®  J  (kp)  cos  (p<|») 

2  n  j-1  J  n-0  p  p 


4  00  S  J  (kp)  sin  (n<)>) 

+  I  v?(  I  ^  °j *  (kal - )  +  Bn  Jn^kP)  sin  Cpd>)  1 

j-1  2  n-1  nUa)  P  P 

n?*p 


where 


< '  rfer  £  1  s!p 


MJij 


BISTATIC  SCATTERING  WIDTH 


The  bistatic  scattering  width  is  called  0  and  is  defined  by 


[15,  eq.  (3-49)] 


0  -  lim  (2ttp|— ||2)  ( 

p-*»  H 

z 

where  H*  is  the  z  component  of  the  incident  magnetic  field  given  by  (1) 

g 

and  Hz  is  the  z  component  of  the  scattered  magnetic  field  given  by 

HS  =  H  -  H1  ( 

z  z  z 

where  is  the  total  magnetic  field  given  by  (56).  Substitution  of 


(56)  into  (70)  gives 


H®  =  H*(M)  +  ESC  -  H*  ,  p  >  a 
z  z  —  z  z 


It  is  evident  from  (1)  and  [14,  eq.  (5-113)]  that 


w 

=  T  e  jnj  (kp)  cos  (n(4>-a)) 
z  A  n  n 
n*0 


Subtracting  (72)  from  (57),  we  obtain 


QA  ,  00  e  j V  (ka)  H^2  (kp)  cos  (n(cf)-a)) 

HSC  -  H1  -  -  l  ~£L— * - TrrS - 

Z  Z  n=0  H(2)  (ka) 

n 


Substitution  of  (60)  and  (73)  into  (71)  gives 


■  J0  [n  (  VjSjn)  -  £n3"Jn(ka)  c°s  (na)  1  ~  H<2) 


H'  J  (kp)  cos  (np) 


+  [  [p  (  [  V°S°  )  -  2jnJ '  (ka)  sin  (not)]  — — pp- 
n-1  n  J-l  j  jn  n  H(2) 


H(2)(kp)  sin  (mt') 


■  *  •  \  >*v 


w;* 

4 

I 


ss 

R 

Jr 

y 


1 


Substituting  (74)  into  (69),  recognizing  that  |H*j  =  1,  using  the 
large  argument  formula  [14,  eq.  (D-13)] 


Hn2)(kp)  * 


and  dividing  both  sides  of  (69)  by  tra,  we  obtain 


£  ■  ±  i  J0  ‘V  <  -  <«>i  ^ 


00  .  n+l  4 


00  .  n+l  4 


+  I  t1—  (  l  V°S°  )  -  2(-l)n  j'(ka)  sin  (net)  ] 
n«l  n  j,1  j  jn  n  H< 


If  there  were  no  slit,  then  (V^  =  V°=*0,  j=l,2,3,4)  and  (76)  would 

SC 

reduce  to  a  quantity  called  a  /ira: 


sc  ,  00  (-1)  e  J'(ka)  cos  (n(<j>-a))  „ 

0 _  m  4  r  _ n  n _  z 

Tfa  Trka  ^ 


H(2)'(ka) 

n 


The  quantity  0  /fra  is  the  normalized  bistatic  scattering  width  of  the 
closed  cylinder,  i.e.,  the  short-circuit  normalized  bistatic  scattering 
width.  This  width  is  that  associated  with  the  well-known  scattering 
pattern  [14,  eq.  (5-117)]  of  the  closed  cylinder. 


VI.  NARROW  SLOT  AT  RESONANCE 


If  the  slot  aperture  is  narrow  and  if  the  associated  cavity  is 
resonant,  that  is,  if  there  is  a  non-negative  integer  p  such  that 
Jp(ka)  *  0,  then  we  can  show  that  the  slit  cylinder  scatters  nearly  the 
same  as  the  corresponding  closed  cylinder  and  we  can  find  approximate 
expressions  for  the  field  components  Inside  the  cylinder.  Thus,  we  will 


have  values  against  which  we  can  check  field  components  and  bistatic 


scattering  widths  computed  from  formulas  in  Sections  IV  and  V. 

In  order  to  do  what  is  proposed  in  the  first  sentence  of  this 
section,  we  consider 


T  Tsc  ,  A 
J  =  J  +  J 


,sc 


where  .J  is  the  electric  current  that  would  be  induced  on  the  cor¬ 
responding  closed  cylinder: 


rsc  -  rusci 

J  =  -  u,  [H  j 

—  — <p  z  p=a 


Moreover,  is  the  electric  current  at  p = a  given  by 
/  »  u.CA® 

-  ~<p  p 


-  u,(Ae  cos  (p<J>)  +  A°  sin  (p<J>)) 


0 


6  O 

where  A  and  A  are  constants.  In  (79)  and  (80),  u,  is  the  unit  vector 
P  P  ~<t> 


in  the  direction.  The  z  component  of  magnetic  field,  and  the  $ 


component  of  electric  field,  E^,  produced  by  radiating  in  the  homoge¬ 


neous  medium  with  constitutive  parameters  y  and  £  are  given  by 


z 


J  (kp)  e  o 

jP(ka)  (Ap  cos  (p<p)  +  Ap  sin  (p$)) 


p  <  a 


(I 


p  >  a 


E^/n 

4> 


j  j;<kp)  e  o 

j  —  (Ap  cos  (p4>)  +  Ap  sin  (p$))  ,  p  £  a 


,  p  >  a 


The  p  component  of  electric  field  produced  by  J_  is  not  considered  here 


Since  E^/0  vanishes  at  p  =  a,  the  $  directed  electric  field 


produced  by  J  of  (78)  cancels  the  $  directed  incident  electric  field 


22 


at  p  ■  a,  Thus,  if  it  did  not  spill  over  into  the  aperture,  J_  would  be 

the  true  electric  current  induced  on  the  slotted  cylindrical  surface. 

e  o  dJd> 

Now,  if  A  and  A  are  adjusted  so  that  both  J,  and  —rf-  vanish  at  <t>  =  0, 

’  P  P  J  4>  dcf) 

then  will  be  small  in  the  aperture  and  therefore  will  be  a  good  approxi¬ 
mation  to  the  true  electric  current.  Here,  is  the  <p  component  of  J_. 
Actually,  if  p  =  0,  then  A°  drops  out  of  £  and  only  can  be  made  to 
vanish  at  =  0.  Substituting  (79)  and  (80)  into  (78),  using  (20),  and 
dJ<t> 

forcing  and  to  vanish  at  $  =  0,  we  obtain 


, ,  00  £  jn  cos  (na) 

I  ” 


TTka 


2H(2)'(ka) 
n=0  n 


(83) 


A°  »  _  _AL. 

p  pirka 


p  ¥  0 


(84) 


In  the  preceding  paragraph,  we  found  that  the  electric  current 

0  o 

on  the  slit  cylinder  is  approximately  J  of  (78)- (80)  with  A  and  A 

-  P  P 

given  by  (83)  and  (84)  if  the  slit  is  narrow  and  if  the  associated  cavity 

is  resonant.  Thus,  the  field  that  exists  in  the  presence  of  the  slit 

cylinder  is  approximately  the  field  which  would  exist  in  the  presence  of 

A 

the  associated  closed  cylinder  plus  the  field  produced  by  J_  .  Since  there 

is  no  field  inside  the  associated  closed  cylinder,  the  field  inside  the 

A 

slit  cylinder  is  approximately  the  field  due  to  J  .  Substituting  (83) 
and  (84)  into  (81)  and  (82)  and  using  [14,  eq.  (D— 17 ) ]  and  the  fact  that 
Jp(ka)  ■  0,  we  express  the  components  (81)  and  (82)  of  this  interior 
field  due  to  as 


00 


H* 

z 


“  £  ]  cos  (na)  1 

2JY'(ka)J  (kp)  [(  l  -P  - )  Cos(p<J>)  +-r  (  l  njnsin(na). 

p  p  n-0  2H(2)  (ka)  p  n-1  „(2)',,  7  sin(p(J,)1 


H(2)  (ka) 
n 


(85) 


>J/n 


00  £  j  cos  (na)  .  00  ,n  ,  .  . 

2Y'(ka)J’(kp)][(  l  - )  cos(pc{>)  +-7  (  [  °-}  *'rW)  sin(p<J>)] 

P  P  ~n  ou(2),^n  P  h(2)  (ka) 


n»0  2H  (ka) 
n 


(86) 


where  the  1/p  terms  are  to  be  discarded  if  p=0.  Here,  Y^(ka)  is  the 
Bessel  function  of  the  second  kind  of  order  p.  The  field  components 

(85)  and  (86)  are  approximately  those  inside  the  slit  cylinder  when  the 

slit  is  narrow  and  the  associated  cavity  is  resonant.  Since,  according 

A 

to  (81)  and  (82),  the  field  due  to  J_  vanishes  outside  the  cylinder,  we 
deduce  that  the  field  outside  this  slit  cylinder  is  approximately  the 
field  that  would  exist  outside  the  corresponding  closed  cylinder. 

In  the  preceding  paragraph,  we  have  shown  that,  if  the  slit  is 
narrow  and  if  the  corresponding  closed  cylinder  is  resonant,  then  the 
field  components  Inside  the  slit  cylinder  are  approximately  (85)  and 

(86)  and  the  field  outside  the  cylinder  is  approximately  the  field  which 
would  exist  if  there  were  no  slit.  Since  neither  (85)  nor  (86)  depend  on 
the  width  of  the  slit,  it  is  evident  that  the  amount  of  TE  field  that  enters 
a  resonant  circular  cylindrical  cavity  through  an  infinitely  long  but 
narrow  slit  does  not  depend  on  the  width  of  the  slit. 


VII.  NUMERICAL  RESULTS  AND  DISCUSSION 

By  means  of  a  computer  program  that  will  be  described  in  a 
forthcoming  report,  |E^/ri|  of  (68)  was  calculated  in  the  aperture,  and 
| E . /H |  of  (68)  and  | H  |  of  (56)  were  calculated  at  the  center  of  the 

y  Z 

aperture  and  at  the  center  of  the  cylinder.  The  normalized  scattering 


widths  a/iTa  of  (76)  and  0  /ira  of  (77)  were  also  calculated  at  (4> , ot)  = 
(0°,  0°),  (180°,  0°),  and  (180°,  180°).  Results  of  these  calculations 
obtained  by  truncating  the  summations  with  respect  to  n  at  n  *  10,000 
are  presented  here. 

The  aperture  field  amplitude  | E^/n [  shown  in  Fig.  4  agrees  well 
with  [1,  Fig.  27]  and  [16,  Fig.  11].  The  data  of  Fig.  4  were  obtained 
by  replacing  (7)  by 

NN 

M®  =  l  V6  M?  ,  NN  <  4  (87 

J-l  J  J 

If  NN  <  4,  the  matrix  Y  and  the  column  vectors  V  and  I  in  (41)  were 

truncated  accordingly.  In  Fig.  4,  the  aperture  field  amplitude  obtained 

by  using  NN  *  1  is  nearly  the  same  as  that  obtained  by  using  NN  =  4. 

However,  the  field  amplitude  obtained  by  using  NN  *  1  is  very  different 

from  that  obtained  by  using  NN  =  4  in  the  wider  slit  aperture  of  Fig.  5. 

Nevertheless,  the  field  amplitude  obtained  by  using  NN  =  3  in  Fig.  5  is 

nearly  the  same  as  that  obtained  by  using  NN  =  4.  The  curve  for  which 

NN  ■  4  in  Fig.  5  agrees  with  [2,  Fig.  4], 

Figure  6  shows  a  ka  scan  of  |E^/n|  at  the  center  (p=a,  4>=0) 

of  the  slit  aperture  in  the  cylinder  of  Fig.  1  for  (p  =  5°  and  a  =  0°. 

o 

The  data  of  Fig.  6  and  of  all  subsequent  figures  were  obtained  by  using 
NN  *  4.  In  Fig.  7,  the  data  of  Fig.  6  are  plotted  on  a  logarithmic 
scale.  The  peak  of  |E^/n|  at  ka  =  0.374  is  probably  due  to  the  H-pol 
low  frequency  resonance  noted  in  [8,  p.  519].  The  three  minima  of 
|E^/ni  occur  very  close  to  resonant  values  of  ka.  Resonant  values  of  ka 
are  values  of  ka  for  which  there  is  a  non-negative  integer  p  such  that 


Jp(ka)  ■  0.  After  dipping  to  a  minimum  very  close  to  a  resonant  value 
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The  solid  curve  in  Fig.  14  is  a  ka  scan  of  the  normalized  back- 
scattering  width  0/TTa  of  the  cylinder  of  Fig.  6.  This  a/TTa  is  given  by 
(76)  with  <J>  ■  a  ■  0°.  The  dashed  curve  in  Fig.  14  is  the  normalized 
back-scattering  width  of  the  corresponding  closed  cylinder,  given  by  (77) 
with  4>  -a*  0°.  Comparing  the  solid  curve  in  Fig.  14  with  that  in  the 
left-hand  plot  of  [8,  Fig.  19b],  we  note  that  the  first  peak  of  the  solid 
curve  in  Fig.  14  is  slightly  higher  than  that  of  the  solid  curve  in  [8]. 
Furthermore,  the  spike  in  the  neighborhood  of  ka  =  5.32  is  more  pro¬ 
nounced  on  the  solid  curve  of  Fig.  14  than  on  the  solid  curve  in  (8]. 

In  Fig.  15,  the  normalized  forward  scattering  width  o/ira  of  the 
cylinder  of  Fig.  6  is  plotted  solid  and  the  normalized  forward  scattering 
width  of  the  corresponding  closed  cylinder  is  plotted  dashed.  These 
scattering  widths  were  calculated  by  setting  4>  =  180°  and  a  =  0°  in  (76) 
and  (77).  Comparing  the  solid  curve  in  Fig.  15  with  that  in  the  right- 
hand  plot  of  [8,  Fig.  19b],  we  observe  that  the  spike  in  the  neighborhood 

of  ka  *  5.32  protrudes  down  farther  in  Fig.  15  than  in  [8], 

The  solid  curve  in  Fig.  16  is  a  ka  scan  of  the  normalized  back- 
scattering  width  a/Tra  of  the  cylinder  of  Fig.  14  with  a  =  180°  instead 
of  a  *  0°.  This  ah ra  is  given  by  (76)  with  $  =  a  =  180°.  The  dashed  curve 
in  Fig.  16  is  the  normalized  back-scattering  width  of  the  corresponding 
closed  cylinder,  given  by  (77)  with  <p  =a=  180°.  In  Fig.  16,  the  aperture 
is  in  the  shadow  region  whereas  in  Figs.  14  and  15,  the  aperture  was 
directly  illuminated.  By  reciprocity  [15,  p.  45],  the  forward  scattering 
width  of  the  cylinder  of  Fig.  16  is  exactly  the  same  as  the  forward  scatter 
ing  width  of  the  cylinder  of  Fig.  14.  Comparing  the  solid  curve  in  Fig.  16 

with  that  in  [8,  Fig.  21a],  we  observe  that  the  first  peak  of  the  solid 

curve  in  Fig.  16  is  slightly  higher  than  that  of  the  solid  curve  in  [8]. 


In  Section  VI,  it  was  shown  that,  at  resonant  values  of  ka, 

scattering  from  a  cylinder  with  a  narrow  slit  aperture  is  approximately 

the  same  as  scattering  from  the  corresponding  closed  cylinder.  This 

means  that  the  solid  and  dashed  curves  in  each  of  Figs.  14,  15,  and  16 

should  cross  each  other  at  values  of  ka  very  close  to  the  resonant  values 

s  c 

of  ka.  Labelling  the  solid  curve  a/ it  a  and  the  dashed  curve  0  /na,  we 

s  c 

computed  0/ Tra  and  a  /TTa  of  Figs.  14,  15,  and  16  at  the  first  six 
resonant  values  of  ka,  given  in  the  first  column  of  Table  1.  Reassuringly, 
the  difference  between  each  of  these  values  of  oh Ta  and  the  corresponding 

SC 

value  of  CJ  /Tra  was  always  less  than  1%.  In  Fig.  14,  o/TTa  is  correct  at 
ka  *  5.31755  and  at  ka  *  5.33144  because  it  agrees  with  o  /rra  there. 
Correctness  of  o/rra  at  these  values  of  ka  in  Fig.  14  is  meaningful  because 
o/t ra  is  changing  so  rapidly  there.  Reading  our  numerical  values  of  a/ Tra 
in  the  neighborhood  of  ka  -  5.32,  we  see  that  a/iTa  is  0.973  at  ka  =  5.31755, 
rises  to  1.49  at  ka  *  5.32667,  dips  to  0.048  at  ka  =  5.326885,  and  then 
rises  to  0.979  at  ka  «  5.33144! 

The  precise  locations  of  the  peaks  of  the  solid  curves  in  Figs.  14, 
15,  and  16  are  pinpointed  in  Table  3.  The  entries  in  any  column  of  Table  3 
are  the  values  of  ka  at  the  peaks  of  the  solid  curve  in  the  figure  whose 
number  appears  at  the  head  of  the  column.  Some  spaces  in  Table  3  were  left 
empty  in  order  to  obtain  alignment  of  the  abscissas  of  similar  maxima. 
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Table  3,  Values  of  ka  at  the  peaks  of  the  solid  curves  in  Figs.  14 
15,  and  16. 
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Fig.  16 

0.375 

0.3725 

0.3755 

0.800 

0.91 

1.958 

1.952 

1.96 

2.30 

2.22 
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3.11 
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3.40 
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4.20 
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APPENDIX  A 


EVALUATION  OF  S?  AND  S? 

Jn  in 


Changing  the  variable  of  integration  from  <f>  to 


x  -  <j>/4> 


(A-l) 


in  (16)  and  (26),  we  obtain 


e  d> 

-e  -  nro 
3jn  "  2n 


,o  !o  ( 

"jn  =  tt  J 


2j-2  /u  . 
x  cos  (bx) 


'1  -  x 


x^  ^sin  (bx) 


'1  -  x' 


(A-2) 


where 


b  =  n4> 


(A- 3) 


(A-4) 


In  regard  to  (A-2)  and  (A-3) ,  it  is  known  that  [17,  formulas  3.752(2.) 


and  3.753(2.)] 


/ - j  (b) 

1  -  x  cos  (bx)  dx  =  — - - 

D 


cos  (bx) 


dx  =  J  (b) 
o 


(A-5) 


(A-6) 


-1  v/1  -  x 


Assuming  that  n  ^  0,  comparison  of  (A-6)  with  (A-2)  gives 


sTn  =  <t>J  (b) 

In  o  o 

Integrating  the  right-hand  side  of  (A-3)  by  parts,  we  obtain 


(A-7) 


4>  b  f 
,o  ro 

5ln  "  j 


'1  -  x  cos  (bx)  dx 


(A-8) 


© 


which,  in  view  of  (A-5),  reduces  to 
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Sln  *  Vl(b) 


(A-9) 


Seeking  a  recurrence  relation  for  when  n  ^  0,  we  use  the  identity 


2  j  —2  2j-4„  ,,  2.. 

X  ■  x  (1  -  (1-x  )) 


to  recast  (A-2)  as 


_e  „e  o 

jn  j-l,n  tt 


2j-4  [7  2  .  . 

x  v/1  -  x  cos  (bx)  dx 


(A-10) 


(A-ll) 


which,  upon  integration  by  parts,  becomes 


S?  =  .  +  —  sin  (bx)  (x2^  ^  J 1  -  x2)  dx 

jn  j-l,n  bit  dx  V 

-1 

Performing  the  indicated  differentiation  in  (A-12) ,  we  obtain 


(A-12) 


S®  "  1  n  ”  "b~31  1  +2(v  ^  S1  2  n*  j=2'3*“*  (A“13) 

jn  j-l,n  b  j-l,n  b  j-2,n 

where  the  S°  0  term  is  to  be  omitted  when  j*>2. 
j  ~2  ,n 

Seeking  a  recurrence  relation  for  S°n>  we  use  an  identity  similar  to 
(A-10)  to  recast  (A-3)  as 


o  „o  o 

jn  “  j-1  ,n  it 


X2J  ^  J 1  -  x2  sin  (bx)  dx 


(A-14) 


An  integration  by  parts  changes  (A-14)  to 


Sjn  =  Sj-l,n  "  ^  \  cos  (bx)  £  (x2j’3  \A  ~  x2)  dx  (A"15) 

-1 

Performing  the  indicated  differentiation  in  (A-15),  we  obtain  the 
recurrence  relation 

SJn  ■  +  ^  SJ„  -  ^  •  J’2’3-"  <A-16> 

Using  (A-13)  and  (A-16)  to  recur  up  from  (A-7)  and  (A-9) ,  we  obtain 
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to? 

wSn 
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W  \rm  1 


bJ  (b) 
,e  _  o 

^ln  n 


e  bJQ(b)  -  JL(b) 


p  (b  -3)  (bJ  -  2J.) 
,e  o  1 


(b4-9b2+60)  bJ  (b)  -  3(b4  -  11  b2  +  40)  J.  (b) 
o  i 


o  bJl(b) 

5ln  n 


(b  -2)  J. (b)  +  bJ  (b) 
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(b4-7b2+24)  J.(b)  +  2b(b2  -  6)  J  (b) 
1  o 


(b6-15b4  +  192b2-720)  J1(b)+3b(b4  -  17b2+120)  JQ(b) 


(A-17) 


(A-18) 


(A-19) 


(A-20) 


(A-21) 


(A-22) 


(A-23) 


(A-24) 


Expressions  (A-19),  (A-20),  (A-22),  (A-23),  and  (A-24)  become  inde- 
ninate  as  b  approaches  zero.  When  b  <  2 ,  we  approximate  cos  (bx)  in 


2)  by  [18,  formula  415.02]. 

2  2  ,4  4  ,6  6  ,8  8  ,10  10 

,  .  bx.bx  bx.  b  x  b  x 

cos  (bx)  =  1  -  2  +  24  ‘  720  +  40320  “  3628800 


ancounter  integrals  defined  by 


(A-25) 


X2j  = 


x2jdx 


-1  \j  1  — x 


,  j  =  0,1,2, , 


is  evident  from  [17,  formula  3.248(3.)]  that 


r  .  "  ■  *  4=n  i  ? 

[2j  (2j ) ! !  *  J  U*  ’  ’ 


.  j  =  0 


(2 j  —1 )  ! ! 


(A-26) 


(A-27) 


(A-28) 


1*3*5  ...  2 j  — 1 ,  j  >  1 


-  j,  .s.'  ’ 
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(2j )  !  !  = 


L  »  j  =  0 

2*  4*  6 ...  2j  ,  j  >  1 


(A-29) 


From  (A-27) ,  we  obtain 
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I2  =  tt/2 


I,  =  3tt/8 
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I.  =  5tt/16 
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Ig  =  35tt/128 
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I1Q  =  63it/256 


I12  -  231tt/1024 
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14 


I16  *  6435tt/32768 


I18  -  12155tr/65536 


In  view  of  (A-30),  substitution  of  (A-25)  into  (A-2)  gives 
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When  b  <  2,  we  approximate  sin  (bx)  in  (A-3)  by  [18,  formula  415.01] 
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